INTERSTELLAR TURBULENCE 
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RESUMEN 

O mcio intcrestelar (ISM) c um sistenia coniplexo, em varias fases, ondc a historia das cstrclas ocorrc. Os 
proccssos dc formagao c morte das estrelas estao fortemente acoplados a dinamica do mcio. Os movimcntos 
caoticos e difusivos obscrvados no meio interestelar caracterizam a turbulencia local. A compreenso da tur- 
bulencia e crucial para o entendimento dos processes de formagao estelar e o feedback de massa e energia ao 
meio a partir das estrelas evoluidas. Campos magneticos tambem permeiam o meio interestelar fazendo dessa 
tarefa ainda mais complicada. Neste trabalho, eu reviso brevemente as principals observagoes e a caracter- 
izagao da turbulencia no meio interestelar. Tambem, mostro os principals resultados teoricos e de simulagoes 
numericas que podem ser usados para construir observaveis, e comparo essas previsoes a algumas observagoes. 

ABSTRACT 

The Interstellar Medium (ISM) is a complex, multi-phase system, where the history of the stars occurs. The 
processes of birth and death of stars are strongly coupled to the dynamics of the ISM. The observed chaotic 
and diffusive motions of the gas characterize its turbulent nature. Understanding turbulence is crucial for 
understanding the star-formation process and the energy-mass feedback from evolved stars. Magnetic fields, 
threading the ISM, are also observed, making this effort even more difficult. In this work, I briefly review 
the main observations and the characterization of turbulence from these observable quantities. Following 
on, I provide a review of the physics of magnetized turbulence. Finally, I will show the main results from 
theoretical and numerical simulations, which can be used to reconstruct observable quantities, and compare 
these predictions to the observations. 
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1. INTRODUCTION 

Turbulence is ubiquitous in the Universe, from 
scales of few centimeters in the water running in the 
sink of our homes, to kiloparsecs at the flowing plas- 
mas in the intergalactic medium (Falceta-Gongalves 
et al. 2010). The classical work of Kolmogorov in 
the 40's have provided a fundamental theoretical ba- 
sis for the study of the dynamics of turbulence. In 
a very simplistic picture, the turbulence is charac- 
terized by chaotic motions in a fluid, resulting in 
diffusion of matter and decay of the kinetic energy 
from large to smaller scales, which end up killed by 
the fluid viscosity with the release of heat. But few 
questions remain: how does the energy cascades from 
the large to the smaller scales? or does this picture 
changes when the fluid is a plasma threaded by mag- 
netic fields? 

In order to answer the first question, Kolmogorov 
assumed that the turbulence in incompressible flu- 
ids should obey the principles of i- homogeneity, ii- 
isotropy, in- scale invariance and iv- locality. By 
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"locality" one should understand wave-wave interac- 
tion at a local wavenumber at Fourier space. In this 
sense, the large scale eddy (fci) decays into smaller 
eddies k2, being k2 = 2ki. The energy transfer 
rate at a given scale / is given by e ~ 5vf/Ti, be- 
ing Ti l/Svi. Therefore, Svi ~ {eiy^^. Regarding 
the power spectrum, if the density is a passive scalar 
we can write Svf E{k')dk\ with k = 27r/Z, 

resulting in the well-known Kolmogorov's spectrum 

In the presence of magnetic fields the turbulence 
cascades quite differently. The main issue is that 
the eddies arc compressed by the field lines in the 
pcrpedicular direction only. In the model presented 
by Goldreich & Sridhar (1995) the timcscales for 
alfvenic and sonic motions are assumed to be similar, 
i.e. ^ Tg, which corresponds to l\\/l± ^ cj^/Svi. 

2/3 

Using the scaling relation of 5vi one obtains /|| oc Ij^ . 
This result shows that the turbulence in magnetized 
plasmas is anisotropic, and that this anisotropy is 
larger at smaller scales. 

Kolmogorov's scalings also fail for compressible 
turbulence. In the case of supersonic fiows the shocks 
are responsible for energy exchanges that results also 
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in compression of the fluid. In this sense, the specific 
energy of a fluid element depends not exclusively on 
its velocity but also on its density. If the energy 
transfer rate is written as e pSv^ / {I / 5v) ~ pSv^/l, 
one may introduce the variable 5u = p^^^6v. The 
power spectrum of Su develops a standard Kol- 
mogorov's slope (i.e. a = —5/3), but a ^ — 2 for 
Sv. 

2. THE TURBULENT ISM 
In the past few decades many important obser- 
vational studies have proven the turbulent nature of 
the interstellar medium (ISM), being i) the distri- 
bution of gas (or column density maps), and ii) the 
dispersion of line profiles, the two major indicators 
of it. 

Armstrong, Rickett & Spangler (1995) presented 
the seminal work in which they used the scintilation 
of background radiation by the turbulent flows in the 
ISM in order to obtain the density spectrum along 
the line of sight. Assuming that the density may be 
considered a passive scalar - it is true for the warm 
and diffuse part of the ISM - the power spctrum may 
also be obtained. Their impressive result was a single 
Kolmogorov slope ranging for more than 6 orders of 
magnitudes in lengthscalc. The data from the den- 
sity distribution in this case corroborated with the 
previous findings of Larson (1981). Larson compiled 
the velocity dispersion distributions observed for sev- 
eral molecular clouds and determined the empirical 
scaling law 6v oc with a ^ 0.3 — 0.5 (on should re- 
member that Kolmogorov's theory predicts a ~ 1/3, 
and for supersonic turbulence a 1/2), as shown in 
Figure 1 (top). 

However, in the past decade new and more sen- 
sitive observations of molecular clouds have shown 
a different picture. When dense cores are plotted 
together with the data compiled by Larson a less 
coherent scaling relation appears. The bottom im- 
age of Figure 1 (extracted from Ballesteros-Paredes 
et al. 2010) shows that the smaller scales present 
a broad range of observed dispersions of velocities, 
and that the Larson relation is not followed by all 
clouds. Ballesteros-Paredes et al. (2010) interpreted 
these observations as an indication that these cores 
were actually contracting. The gravitational collapse 
would result in an increase of the velocity dispersion 
at the scales of the denser cores. 

Magnetic fields can also be used in order to study 
the ISM turbulence. Polarization maps - in the op- 
tical for dust absorption, and in the infrared by dust 
emission - are useful in unreaveling the topology of 
the magnetic field in the plane of sky. The per- 
turbations in the general trend of the field lines. 
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Fig. 1. Dispersion of velocity vs lengthscale for typical 
large molecular clouds (top) - original Larson relation - 
and for all structures, including the dense massive cores 
(bottom) (extracted from Ballesteros-Paredes 2010). 

i.e. the uniform component, are usually associated 
to the turbulent motions in the plasma. In the 
Chandrasekhar-Fermi approximation the equiparti- 
tion between the magnetic perturbation and the tur- 
bulent kinetic energy is assumed. In this case it is 
possible to estimate either the magnetic pressure or 
the dispersion of velocity once the other is known. 

2.1. Results from numerical simulations 

One of the main difficulties faced by observers 
in studying turbulence from the observations is its 
three dimensional nature. Observations arc either 
1-D, when projections along the line of sight are 
considered (e.g. dispersion of velocities, scintilation 
of background radiation) , or 2-D when maps in the 
plane of the sky are used (e.g. column density statis- 
tics and polarization maps). The conversion from 
the observed quantities to the actual physics of tur- 
bulence is not straightforward (see Burkhart et al. 
2009). In this sense numerical simulations represent 
a powerful tool. It is possible to simulate three di- 
mensional turbulent flows and reconstruct obscrv- 
ables from these cubes, considering any position for 
the observer. 

In Falceta-Gongalves, Lazarian & Houde (2010) 
we calculated synthetic observations of line profiles 
and presented dispersion of velocity vs lengthscale 
relationship for different turbulent regimes (Figure 
2), being subsonic (left) and supersonic (right). It is 
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Fig. 2. Dispersion of velocity vs lengthscale, obtained for all LOS's with beam sizes / x I, for two 3D-MHD numerical 
simulations. Left: subsonic/subalfvenic and right: supersonic/subalfvenic (extracted from Falceta-Gongalves, Lazarian 
& Houde 2010). The line corresponds to the actual three dimensional scaling, i.e. l^. 



clear in these calculations, even more in the super- 
sonic case, that the range of velocity dispersions at 
small scales is very large. Is it gravity the main cause 
for that? Actually, these simulations were MHD, 
isothermal, non selfgravitating models. The main is- 
sue is that the "small scale" in the plot does not truly 
corresponds to the smaller turbulent cells. Since the 
line profiles are obtained integrating along the line 
of sight, it is obtained in a volume V — I x I x L, be- 
ing L the largest scale. Within the LOS, many small 
dense structures - each being subject to a drift veloc- 
ity of the order of the largest scale amplitude - are 
intercepted and contribute to the total dispersion. 

The origin of the condensations is related to 
the turbulence sonic Mach number. It reveals that 
shocks are the most important mechanism in gen- 
erating dense structures in the ISM (Ballesteros- 
Paredes, Vazquez-Semadeni & Scalo 1999). 

Regarding the magnetic fields, numerical simula- 
tions have shown that even subalfvenic turbulence, 
at large scales, may present superalfvcnic flows at 
smaller scales (Falccta-Goncalvcs, Lazarian & Kowal 
2008). When shocks occur predominantly perpen- 
dicular to the mean field lines the density increases 
but the magnetic pressure not. Within the denser 
regions, the alfven speed (as it scales with p~^^^) 
decreases and may become smaller than dvi . The su- 
peralfvcnic turbulence at the small scales is respon- 
sible for the bending of the field lines and, e.g. the 
decrease in the polarization degree, as observed in 
dense molecular cores. 



3. IMPLICATIONS ON STAR FORMATION 

It is well stablished that the ISM turbulence plays 
a key role in the star formation process, together 
with the magnetic field. In the same time, stellar 
energy feedback (mostly SNc) triggers the turbulence 
at large scales. Therefore, a full understanding of 
the turbulence cascade and the evolution of magnetic 
fields in turbulent flows is mandatory. 

Uniform magnetic fleld lines are able to prevent 
the gravitational collapse of a cloud in the perpen- 
dicular direction, though not efficient in preventing 
the collapse as a whole. Alfvcnic motions within the 
collapsing clouds may be effective in slowing down 
the inward fluxes of gas and, therefore, reduce the 
star formation efficiency (Falceta-Goncalves, de Juli 
& Jatcnco-Pcrcira 2003). 

Another important issue is the survival of the 
dense cores to the turbulence itself. In one hand, 
the turbulent motions result in shocks and the for- 
mation of the dense clouds. In the other hand, the 
turbulence is responsible for the diffusion of the gas 
and the disruption of cloudlets and cores. How long 
can a star forming region survive? 

From the fragmentation theory for star forma- 
tion, a large molecular cloud gets gravitationally un- 
stable and contracts. During the contraction its den- 
sity gets larger and larger resulting in a decrease 
of the Jeans mass. Parts of this cloud with masses 
larger than the new critical mass "dettach" from the 
whole cloud and starts to collapse by its own. The 
process occurs until the optical depth of the cloud be- 
comes much larger than unity. The last fragmented 
cores represent the final stage of a prestellar cloud. 
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Fig. 3. Clump lifetimes given in terms of the sonic cross- 
ing time. Dashed line corresponds to tufc — tcroas and 
dotted line corresponds to tufo — lOtcross. (extracted 
from Falceta-Gongalves & Lazarian 2011, submitted) 

The fragmentation and collapse of a cloud is, in gen- 
eral, believed to be fast. Much faster than other 
typical timescales of the system. Indeed, the ratio 
between the free-fall time and the turbulent turnover 
timcscale for a Giant Molecular Cloud (GMC) is: 
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However, a typical Ipc sized dense core within 
a molecular cloud presents tff /tturb ~ 10, i.e. the 
turbulence is able to fully develop during the core 
collapse. In order to test the survival of clumps 
to the turbulent diffusion we performed a standard, 
512'^ resolution, 3D MHD numerical simulation with 
Ms ^ 7, and identified and followed clumps for a 
long time. In Figure 3 we present the core lifetimes 
vs the sonic crossing time for all identified overden- 
sity structures (clumps) in the run. It is clear that 
most of the detected clumps survive more than one 
crossing time, but only few can survive more than 
10 crossing times (similar to the values found by 
Galvan-Madrid et al. 2007). The turbulent diffu- 
sion is able to disrupt a clump before it fragments 
and collapses even further. Less than 5% of the 
clumps presented lifetimes larger than 10 crossing 
times. This result presents a possible explanation 
for the decreased star formation efficiency, even in 
regions where there is no massive star or any other 
source of strong radiation. 

A rough estimate of the clump lifetimes is ob- 
tained assuming that all the gas shocking on a clump, 
with size I and density pi, is trapped: 



where z;l is the large scale velocity and p is the den- 
sity of the unperturbed gas. This results in: 
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Equation 3 shows that the lifetime of a small scale 
clump depends on the large scale turbulent eddies as 
well. Also, provided that the column density may be 
obtained by E ~ p;^, we get the scaling tufc oc E. 
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